Higgs boson decays to 77 and Z7 in models with Higgs extensions 
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The decays of a Higgs boson to the 77 and Z'y final states are purely quantum mechanical 
phenomena that are closely related to each other. We study the effects of an extended Higgs sector 
on the decay rates of the two modes. We propose that a simultaneous determination of them and the 
ZZ mode is a useful way to see whether the Higgs boson recently observed by the LHC experiments 
is of the standard model type or could be a member of a larger Higgs sector. 
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Introduction - The quest for the origin of elementary 
particle masses is arguably one of the most important 
tasks in current high energy physics. According to the 
standard model (SM), a scalar field is employed to break 
the electroweak (EW) symmetry down to the electro- 
magnetic (EM) symmetry, SU{2) L x U{l) Y -> C/(1) E m, 
giving masses to the W and Z bosons that mediate 
weak interactions. This so-called Higgs mechanism [l| 
is achieved when the scalar field spontaneously acquires 
a nonzero vacuum expectation value (VEV) due to the 
instability in its potential. As a consequence, the SM 
predicts the existence of a spin-0 Higgs boson. With the 
introduction of Yukawa interactions, fermionic particles 
can obtain their masses from the same Higgs field as well. 
Therefore, the discovery of the Higgs boson does not only 
complete the particle spectrum of the SM, but also re- 
veals the secrets of EW symmetry breaking and mass. 

Recently, both ATLAS and CMS Collaborations @, 0] 
of the CERN Large Hadron Collider (LHC) have reported 
the observation of a Higgs boson at around 125 GeV 
at ~ 5cr level through the combination of ZZ and 77 
channels. Minor excesses are observed in the WW and 
Wh/Zh with h — > bb channels, while no excess is seen in 
the t + t~ channel. Around this mass region, another de- 
cay channel that is closely related to the diphoton mode 
and clean in the LHC environment is the Zj mode. At 
the leading order in the SM, both the 77 [3] and Z7 [f| 
are loop processes mediated by the same particles. New 
particles beyond the SM can change their relative magni- 
tudes. Although the ZZ decay occurs at tree level and is 
less sensitive to new particle contributions, the rate de- 
pends on how the EW symmetry is broken. Therefore, a 
simultaneous measurement of their production rates will 
be helpful in diagnosing the observed Higgs boson. 

In view of the 125 GeV Higgs boson, denoted by h, we 
consider models that have only an extended Higgs sec- 
tor for simplicity. There are some recent studies in the 
literature about the h — > 77 decay in models with an ex- 
tended Higgs sector In this Letter, we investigate 
both the h — > 77 and Z7 decays in models with Higgs 
extensions based on various physics motivations. We first 
assume that h is SM-like, meaning that the couplings of 
h with fermions (hff) as well as the weak gauge bosons 
(hVV) are the same as the SM ones. This situation can 



be realized by imposing a discrete Z2 symmetry, where 
the SM particles are ^-even while the extra scalar fields 
are i^-odd. In this case, the dominant production mech- 
anism of h at the LHC is the gluon fusion process, as in 
the SM. We study how the decay rates of h — > 77 and 
Z7 are modified (see also Q). 

Since no excess is observed in the h — > t + t~ mode and 
the Wh/Zh, h — > bb decay has only a minor excess, it is 
still unclear whether the observed Higgs boson is SM-likc. 
We therefore also analyze the scenario of a fermiophobic 
(FP) Higgs boson, in which the Yukawa couplings are 
vanishingly small. In this case, h is dominantly produced 
via the weak boson fusion process [l(| ■ 

Models with a Higgs extension - In general, models 
with an extended Higgs sector often contain charged 
Higgs bosons that, among other phenomena, can con- 
tribute to the h — > 77 and Z7 decays through the loop 
effect. Although there are many possibilities for the ex- 
tended Higgs sector, we discuss those with extra SU(2)l 
singlets S (with Y = 1 and Y = 2), doublet D (with 
Y = 1/2), and triplets T (with Y = and Y = 1), whose 
charge assignments are given in Table H] [lTj . 

We consider three distinct classes of extended Higgs 
sectors: (I) models with one singly-charged scalar bo- 
son, (II) those with one singly-charged and one doubly- 
charged scalar bosons, and (III) those with two singly- 
charged scalar bosons. According to the representations 
listed in Tablc[IJ there are three, four, and six possibilities 
for Class I, Class II and Class III, respectively, all listed 
in Table [TlJ Examples of models in Class I (Models 1 to 
3) include the two Higgs doublet model [ll[ and the min- 
imal supersymmctric SM. Models in Class II (Models 4 
to 7) include the Higgs triplet model [l2[ and Zee-Babu 
model [I3|. Finally, models in Class III (Models 8 to 13) 
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TABLE I: Charge assignments for extra scalar fields under 
the SU(2)l x U(1)y gauge symmetry. 
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TABLE II: Models considered in this work and their extra 
scalar field contents. 



sw = sin 9w, cw — cosOyy with 9w being the weak mix- 
■ ing angle. In Class II models, Mj for the singly-charged 
scalar boson (i = 1) and for the doubly-charged scalar 
boson (i = 2) are generally independent parameters, ex- 
cept for Model 5 where they are the same. 

In Class III models, on the other hand, the two singly- 
charged charged scalar bosons Sf and Sf generally mix 
with each other, so that the expressions for g l g SZ and 
" X l g Sh can be quite different from those given in Eq. (|3j). 
In this case, the coupling g l g SZ are written in the mass 
cigenbasis of the two charged scalar bosons (Sf, S 2 ) as 



include those where tiny Majorana neutrino masses are 
generated via higher- loop processes [14 1 . 

Loop-induced diboson decays of Higgs - The decay 
rates of h — > 77 and Zj can be generally expressed as 
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where Gf is the Fermi decay constant, v — l/(v / 2Gi?) 1 / 2 
is the Higgs VEV, rrih is the Higgs boson mass, mz is the 
Z boson mass, Qx is the electric charge of particle X, N[ 
is the color factor of the fermion /. The loop functions 
for the scalar contribution P s and Jg are given by 

I l s = — [1 + 2m s .C (0,0, to 2 , m Si ,m s% ,m Si j\ , 
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in terms of the Passarino-Veltman functions Bq and Co 
defined in [l5[, where mgj is the mass of the charged 
scalar boson Si- The loop functions for the W boson 
(Iw and Jw) as well as the fermion / (If and Jj) con- 
tributions are given in [l6| . The couplings between the 
charged scalar bosons Si and the Z boson as well as h 
are derived from 



-\% Sh SiS*h + ig % i sz (d^SiS* + Sid^ZV + h.c. 
In models of Classes I and II, the couplings 



9ssz 



f(4 

Cw 



SwQst), X^ sh = -(mi-Mf) , (3) 



where M, is a dimensionful parameter unrelated to the 
Higgs VEV, I\ is the third isospin component of Si, and 
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where 9 is the mixing angle (eg = cos 9, sg = sin 6*) con- 
necting between the weak eigenstates (tpf,tp 2 ) and the 



mass eigenstates: 
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^ i^ 2 , the off-diagonal couplings g Z ggz are 
nonzero and can contribute to the h — > Zj decay as well. 
A< SSh can ^ e calculated for Models 8, 9 and 10 as 
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Those for Models 11 and 12 are 
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with i = 1,2 . (5) 



In the above expressions for Ag S/l , the dimensionful pa- 
rameters Mi 2,3 show up in the scalar potential 
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where ip\ and (^2 are the scalar fields including ipf 1 and 
ip 2 , respectively. In Models 11 and 12, the parameter cor- 
responding to M3 is absent, while there is another dimen- 
sionful parameter fi defined in the terms fi<f>' cpi<p2 + h.c. 
that induce mixing between iff and <p 2 , where <!> is the 
Higgs doublet field associated with h. In Model 13, there 
are no parameters corresponding to M3 and y, and, there- 
fore, there is no mixing at tree level. 

We note in passing that the coupling formulae for \ l g Sh 
are derived assuming an unbroken discrete Z 2 symmetry 
in each model. If the Z 2 symmetry is broken or does not 
exist at all, Ag Sfl are generally different from those given 
above, even when we take the limit where h is SM-likc. 

Numerical results - To see the correlation between the 
decay rates of h — > Z7 and 77, we further define the 
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Classes I, II and III, respectively: 



TABLE III: Total decay rate and branching ratios of a regular 
and FP Higgs boson, respectively, in the SM. 



Model 


1 


2 3 




A23 77 (%) -9.9 (4.2) -9.9 (4.2) -9.9 (4.2) 
ABz-y (%) 3.3 (-1.2) -3.8 (1.4) -11 (4.2) 
R 0.72 (0.60) 0.67 (0.61) 0.63 (0.63) 


Model 


4 


5 6 
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AZ3 77 (%) 
AB Zl (%) 
R 


-45 (11) 
17 (-3.4) 
1.3 (0.55) 


-45 (11) -45 (11) -45 (11) 
-12 (2.5) 9.4 (-1.9) 2.2 (-0.50) 
1.0 (0.58) 1.3 (0.56) 1.2 (0.56) 



TABLE IV: Deviations of the h — > 77 and Z7 branching 
ratios from the SM expectations and the ratio R in models 
of Classes I and II. All the charged scalar boson masses are 
taken to be 400 GeV. The numbers without (with) paren- 
theses correspond to the case with M — (M — Mmax), 
where Mmax denotes the maximum value allowed by the con- 
ditions in Eqs. © and (|7|. In Class I models Mmax = 470 
GeV, while in Class II models Mmax = 440 GeV assuming 
Mi = M 2 = M. 



Model (FP) 


1 2 3 




AZ3 77 (%) 
A23z 7 (%) 
R 


500 (570) 500 (570) 500 (570) 
360 (340) 330 (350) 300 (360) 
0.48 (0.41) 0.45 (0.42) 0.42 (0.43) 




Model (FP) 


4 5 6 


7 


A23 77 (%) 
ABz-, (%) 
R 


320 (610) 320 (610) 320 (610) 320 (610) 
410 (330) 290 (350) 380 (330) 350 (340) 
0.78 (0.38) 0.60 (0.40) 0.73 (0.39) 0.69 (0.39) 



TABLE V: Same as Table [TV] but for the FP Higgs case. 

ratio R = r^ 7 /r 77 . First, we give the SM expectations 
of the two diboson decays of a regular and FP Higgs 
boson in Table IIII| where = 125 GeV and m t = 
173 GeV are used. Next, we show numerical results for 
the case with an extended Higgs sector. As mentioned 
before, we assume that the observed Higgs boson is SM- 
like in couplings with the gauge bosons and fermions in 
our numerical studies. Moreover, we present the results 
for the case where the charged scalar boson are at least 
400 GeV in mass. 

For meaningful discussions, we calculate parameter 
bounds by considering unitarity and vacuum stability 
constraints. For unitarity, we impose the condition that 
the magnitudes of all dimensionless coupling constants do 
not exceed 2n. For vacuum stability, we require that the 
scalar potential is bounded from below in the parame- 
ter space where the quartic terms dominate. Combining 
the two bounds, we obtain the following conditions for 
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We shall use them and neglect the slight changes among 
the models in each class due to some coupling constants 
that are irrelevant to our analysis. 

The deviations from the SM predictions for the h — >• 77 
and Z"f branching ratios can be parametrized as 
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where B^ (Bz~) is the branching ratio of h — > 77 (h—> 
Zj) in a Higgs-extended model, while B^f (Bff) is that 
in the SM. In Table IIV( the deviations and R are listed 
in the case of m s + = 400 GeV for Class I models, and 
m s + = m s ++ = 400 GeV with M 1 = M 2 = M for 
Class II models. The parameter M is constrained by 
Eqs. © and to be < M < 470 (440) GeV for Class I 
(Class II). The results for the FP Higgs boson scenario 
using the same parameter set are given in Table |VJ 

For models in Class III, we consider as an example 
the case where m s + = 400 GeV and m s + = 500 GeV. 
We assume that the three dimensionful parameters are 
the same: Mi = M2 = M3 = M. The maximally al- 
lowed value of M depends on the mixing angle 6, but the 
strictest upper bound on M from Eq. ((Sj) is found when 
9 = 0. In this case, the upper bound is 440 GeV, which 
is used in the following numerical analysis. 

In Fig. [H the deviation of the h — > 77 branching ratio 
is plotted against sin£>. According to the left plot for 
the scenario with a regular Higgs boson, Models 8-10 
all have the same behavior of positive corrections, while 
Models 11 and 12 share a common negative corrections. 
Model 13, on the other hand, is a constant due to Eq. ([5]) 
and has virtually no deviation from the SM expectation. 
In the scenario with an FP Higgs, however, the deviations 
arc all shifted to around +550% or more. 
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FIG. 1: Deviation of the branching ratio of h — > 77 as a 
function of sintf. We take Mi = M 2 = M 3 = M = 440 GeV, 
m q + = 400 GeV and m q+ = 500 GeV. The left (right) figure 
shows the case with the regular (FP) Higgs boson. 
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FIG. 2: Deviation of the branching ratio of h — ► Z7 as a 
function of sin 6*. We take Mi = M 2 = M 3 = M = 440 GeV, 
m„+ = 400 GeV and m„+ = 500 GeV. The left (right) figure 
shows the case with the regular (FP) Higgs boson. 



= 50(1 GeV. M = 440 GeV 



511(1 GeV. M = 440 GeV 




FIG. 3: The ratio R as a function of sin 9. We take Mi = 
M 2 = M 3 = M = 440 GeV, m„+ = 400 GeV and m„+ = 500 

GeV. The left (right) figure shows the case with the regular 
(FP) Higgs boson. The SM case is indicated by the dashed 
line for comparison. 

The deviation of the h — > Z7 branching ratio varies a 
lot among the Class III models, as shown in Fig. [3J For 
the scenario with a regular (FP) Higgs boson, the devi- 



ation is between about —3% and 5% (330% and 360%). 
Finally, we show the ratio R in Fig. [3] for the two scenar- 
ios as well, in contrast with the values for the SM. 

We now briefly comment on the case with lighter 
charged scalar bosons. When we take all their masses 
to be 200 GeV and use the maximally allowed M, the 
deviation in the branching ratio of h — > 77 is predicted 
to be 14 %, 59% and 22% in Classes I, II and III with 
9 = 0, respectively, in the regular Higgs boson scenario. 

Summary - With the observation of a Higgs boson h 
of mass 125 GeV by ATLAS and CMS, it would be in- 
teresting to diagnose whether h is standard model-like 
(SM-like) or part of a larger Higgs sector. We thus con- 
sider and classify models with simple Higgs extensions. 
We have imposed the unitarity and vacuum stability con- 
straints on the model parameters. We have studied the 
neutral diboson decays of h, assuming that it has SM- 
like couplings with the weak bosons and fermions. In our 
framework, the ZZ mode is virtually unaffected, whereas 
the 77 and Zj modes can be modified by a few to a few 
tens of percent. A simultaneous determination of their 
branching ratios is thus useful in exploring the possibil- 
ity of an extended Higgs sector. As a comparison, we 
have also considered the scenario where h is fcrmiopho- 
bic and observed more dramatic changes in the diboson 
decay rates. 
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